Microwave spectroscopy of a Cooper pair beam splitter 



Audrey Cottet 

Lahoratoire Pierre Aigrain, Ecole Normale Superieure, CNRS UMR 8551, 
Laboratoire associe aux universites Pierre et Marie Curie et Denis Diderot, 
24, rue Lhomond, 75231 Paris Cedex 05, France 
(Dated: March 3, 2013) 

This article discusses how to demonstrate the entanglement of the split Cooper pairs produced 
in a double-quantum-dot based Cooper pair beam splitter (CPS), by performing the microwave 
spectroscopy of the CPS. More precisely, one can study the DC current response of such a CPS 
to two on-phase microwave gate irradiations applied to the two CPS dots. Some of the current 
peaks caused by the microwaves show a strongly nonmonotonic variation with the amplitude of the 
irradiation applied individually to one dot. This effect is directly due to a subradiance property 
caused by the coherence of the split pairs. Using realistic parameters, one finds that this effect has 
a measurable amplitude. 

PACS numbers: 73.23.-b,73.63.Fg,03.67.Bg 



I. INTRODUCTION 

Quantum entanglement between spatially separated 
particles represents a promising resource in the field of 
quantum computation and communication. However, 
this fascinating behavior can be difficult to observe in 
practice due to decoherence caused by the particles en- 
vironment. This is why the "spooky action at a dis- 
tance" was first demonstrated with photons, atoms, or 
ions which can be naturally placed in weakly interacting 
conditions l|- 



Observing electronic entanglement in solid state sys- 
tems is a-priori more challenging since an electronic fluid 
is characterized by a complex many-body state in gen- 
eral. However, quantum entanglement has been recently 
observed on superconducting chips [4]. In this case, the 
particles are replaced by superconducting quantum bits, 
which can be sufficiently well isolated from the out- 
side world thanks to the rigidity of the superconducting 
phase, if an appropriate circuit design is used. In these 
experiments, the entangled degrees of freedom are de- 
fined from the charges of small superconducting islands, 
or from the persistent current states of a superconducting 
loop, for instance [5]. 

Superconductors enclose another natural source of en- 
tanglement which has not been exploited so far, i.e. the 
spin entanglement of its Cooper pairs. In a conventional 
superconductor. Cooper pairs gather two electrons cor- 
related in a spin-singlet state. The use of this resource 
for entanglement production requires to build hybrid cir- 
cuits in which the superconductors are connected to non- 
superconducting elements which allow the spatial separa- 
tion of Cooper pairs. In principle, a double quantum dot 
circuit connected to a central superconducting contact 
(input) and two outer normal metal contacts (outputs) 
facilitates this process [6]. Such a "Cooper pair splitter" 
(CPS) has been realized recently by using double dots 
formed inside semiconducting nanowires [7-9] or carbon 



nanotubesfiol. |Tl|. The spatial splitting of the Cooper 
pairs has been demonstrated from an analysis of the cur- 
rent response of the CPS to a DC voltage bias. However, 
the spin entanglement of the split pairs was not tested 
by these experiments. 

It has been suggested to use the noise cross-correlations 
of the electrical current to characterize the degree of en- 
tanglement of pairs of electrons |12l - l2l| . Alternatively, 
Ref. (22! proposes to put in evidence spin entanglement 
by coupling the CPS to a microwave cavity. In this ref- 
erence, a double quantum dot formed inside a single wall 
carbon nanotube is considered. Spin-orbit interaction 
produces a coupling between electronic spins and cavity 
photons. Such a coupling leads to a lasing effect which in- 
volves a transition between the spin singlet state in which 
Cooper pairs are injected and some spin triplet states. 
This effect vanishes when the spin/photon coupling is 
equal in the two dots, due to a subradiance property 
caused by the entangled structure of the spin-singlets. 
However, realizing such an experimental scheme is chal- 
lenging since it requires to couple a complex quantum 



dot circuit to a photonic cavitv|23l-l2 



The present work suggests an alternative strategy to 
exploit the subradiance of spin-orbit induced transitions 
between spin singlet and spin triplet CPS states. One can 
measure the DC current at the input of the CPS when 
microwave gate voltage excitations are applied separately 
to the two CPS dots. The microwave-induced state tran- 
sitions mediated by spin-orbit coupling result in current 
peaks at the input of the CPS versus the dots DC gate 
voltages. Assuming that two on-phase microwave exci- 
tations are applied to the two dots, these peaks vanish 
when the amplitude of the two excitations become equal. 
This subradiant behavior is directly related to the spin- 
entanglement of the split Cooper pairs hosted by the 
CPS. 

This article is organized as follows. Section II defines 
the CPS hamiltonian, for a single wall carbon nanotube 
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based implementation. Section III discusses the CPS 
even-charged eigenstates in the absence of the microwave 
excitations and without the normal metal contacts. Sec- 
tion IV describes the coupling between the CPS even- 
charged eigenstates and the microwave excitations. Sec- 
tion V describes the CPS state dynamics in the presence 
of the voltage-biased normal metal contacts, by using a 
master equation description. Section VI describes the 
results given by this approach, an in particular the pre- 
dictions obtained for the DC current at the input of the 
CPS. Section VII presents further examination and mod- 
ifications of the model, which are useful to put the results 
of section VI into perspective. In particular, it discusses 
the role of atomic-scale disorder in the nanotube, the role 
of the form assumed for the spin-orbit interaction term, 
and possible microwave induced transitions in the CPS 
singly occupied charge sector. Section VIII compares the 
measurement strategy discussed in this work to the one 
of Ref. [22|. Section IX concludes. Although this article 
focuses on a carbon-nanotube-based CPS, the entangle- 
ment detection scheme discussed in this work could be 
generalized to other types of quantum dots with spin- 
orbit coupling like e.g. InAs quantum dots, in principle. 
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FIG. 1: Scheme of a Cooper pair splitter made out of a carbon 
nanotube. The two quantum dots L and R are defined by the 
normal metal contacts (in green) and the superconducting 
contact (in blue) deposited on top of the carbon nanotube 
(in light blue). The dot L{R) is capacitively coupled to a 
DC gate voltage V^^^^ and microwave gate voltage Va}^^ - 
The superconducting contact is connected to ground and the 
normal metal contacts are biased with a voltage Vh. 



II. HAMILTONIAN OF THE CPS 

Let us consider the circuit represented schematically in 
Fig. 1. Two normal metal contacts and a superconduct- 
ing contact are used to define two quantum dots L and R 
along a single wall carbon nanotube. The superconduct- 
ing contact is connected to ground, and a bias voltage 
Vh is applied to the two normal metal contacts. The 
dot L{R) is connected capacitively to a DC gate volt- 
age source V^^^^ and a microwave gate voltage source 
Vac^\t). In the following, it is assumed that V^^{t) 
and V^(t) are in-phase, i.e. Vac^\t) = Vac^^ sin(cji^i?t). 
Inside the left and right dots i G {L, i?}, an electron with 
spin cr G {t^i} can be in the orbital r G {K^K'} of the 
nanotube, which is reminiscent from the K / K' degener- 
acy of graphene. One can use a double dot hamiltonian 
which takes into account the proximity effect caused by 
the superconducting contact, i.e. 



^DQD = X] + /^sor(7)nira + Hprox (1) 



with 



orbital r of dot i 



d\^^ the creation operator for an electron with spin a in 



dl^^dira- For 



G {L, R} and riircr 
simplicity, one can assume that the orbital energies in 
dots L and R are both equal to £ in the absence of the 
external microwave irradiation, which can be obtained 

Li IV) 

by tuning properly the dots' DC gate voltages Vg . 
The term A^o is caused by spin-orbit coupling inside the 
carbon nanotube [26|. The term /S^k^k' describes a cou- 
pling between the K and K' orbitals of dot i, due to 
disorder at the level of the nanotube atomic structure 



29] . The term in tee describes interdot hopping. The 
term Hint accounts for Coulomb charging effects. One 
can assume that there cannot be more than one electron 
in each dot, due to a strong intra-dot Coulomb charging 
energy. Therefore, Cooper pairs injected inside the CPS 
are split into two electrons, one in each dot. The term 
Hprox accounts for coherent injection of singlet Cooper 
pairs inside the double dot [30]. This approach is valid 
provided quasiparticle transport between the supercon- 
ducting contact and the double dot can be disregarded. 
This requires eVb < A, with A the BCS gap of the su- 
perconducting contact. The hamiltonian Hdqd must be 
supplemented by the normal leads hamiltonian 



leads 



(3) 



and the tunnel coupling between the dots and normal 
leads 



h.c. 



(4) 



with Cik^a the annihilation operator for an electron with 
spin cr in orbital kr of the normal lead i G {L, R}. 
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The effect of the microwave gate voltage bias can also 
be described with hamiltonian terms. The gate voltage 
Vac^\t) = Vac^^ ^ni{ujRF't) corrcsponds to an electric 
field Eac^^ = Vac^\t) / with d the center to ground 
separation of the waveguide providing the microwave sig- 
nal. This also corresponds in the Coulomb gauge to a 

L(R) L(R) 

vector potential Aac — —Vac cos{ujRFt)/ujRFd on dot 
L{R), which is assumed to be perpendicular to the carbon 

L/( R) 

nanotube. The interplay between Aac and inter sub- 
band spin-orbit coupling elements induced by the nan- 
otube curvature results in a spin/photon coupling term 
(see Ref. [31] for details) 



^RF 



■E 



, cos{u RFt)dl^^d 



(5) 



with e > the electron charge. For simplicity, this arti- 
cle uses the particular structure aira = '^crai with G M 
and i the imaginary unit number, obtained from a micro- 
scopic description of spin-orbit coupling in a zigzag nan- 
otube quantum dot [3l|, based on Refs.^, ^] (see also 
Refs. I34l438l ). However, part VII. 1 will show that the re- 
sults presented here can be generahzed straightforwardly 
to a more general a^ra- The dimensionless coefficient 
corresponds to the coefficient Xi/eVrms of reference [22|, 
with Vrms the amplitude of vacuum voltage ffuctuations 
for the photonic cavity considered in this reference. The 
value of ai can be estimated to typically 3.10"^ while 
Vac^^ can reach typically 100/iV. One can also use a 
hamiltonian term Hf^p to account for the modulation of 
the dots orbital energies by the microwave gate voltages. 
For simplicity, one can disregard the mutual capacitive 
coupling between the two dots. In this case, one finds 



^RF 



-^Kicvl^ sin{ujRFt)ni. 



(6) 



where hCi is a dimensionless capacitive coupling constant 
which is typically of the order of 10~^. 

In the following, it is assumed that electrons can go 
from dot L{R) to the corresponding normal metal contact 
but not the reverse. This can be obtained by using a bias 
voltage Vb such that 



eVb > 2Ar + te 



{S-2Ar)^^XkBT (7) 



with Ar = \J A^^ + ^K^K' ^ ^ dimensionless co- 
efficient which takes into account the effective thermal 
broadening of the levels (see Ref. [31] for details). 



^ = 26 the energy of a CPS doubly occupied state for 
te/i = Aso = Ak^k' — 0. The parameter ^ can be tuned 

Ij{FO\ 

with Vg • The coupling tg/i hybridizes the CPS empty 
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FIG. 2: (a) Energies Ex, Ei and E- of the states l^i), \V2) 
and |T_) as a function of b. (b) Transition frequencies uv-^t_ , 
l>jt_V'2 of the CPS as a function of 8. (c) Dynamics 

of the CPS near the working point 8 — 2Ar. We consider a 
bias voltage regime such that the tunnel transitions between 
the different CPS states (blue arrows) occur together with the 
transfer of one electron towards the normal contacts. A mi- 
crowave irradiation can induce transitions between the states 
\Vx) and \V<i), \Vi) and |T_), or |T_) and 1^2) without any 
transfer of electrons between the CPS and the leads (red wavy 
arrows). We have used teh/Aso = 1/3 and Ak/k' / Aso = 6 
in panels (a) and (b). 

state |0, 0) with the subspace of the CPS doubly occupied 
states {|rcr, rV)}, where \tg^t'cf) denotes a CPS state 
with one electron with spin cr in orbital r of dot L and 
one electron with spin a' in orbital t' of dot R 
The resulting even-charged subspace is called £. Near 
the working point 5 ^ 2Ar, the CPS dynamics involves a 
subspace S of at maximum five eigenstates from S. Three 
of these eigenstates have an energy E- = J — 2A^, namely 



(8) 



2Ar 



^|C+(r t,r;)) 



III. EXPRESSION OF THE EVEN-CHARGED ,^ . ^ 1 / A^^ \ \K(j,K(j) - \K'(j,K'(j) , , 



CPS EIGENSTATES 1" + / ~ 2 V A^ / 

This section discusses the relevant eigenstates of +^^(7 \Cj^{Kcr^K' cr)) 

-^DQD even charge sector for 5 ~ 2Ar, with o- ^ 
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and 

(10) 

where a denotes the spin direction opposite to a and 
|C±(rcr,rVO) = (|rcr,rVO ± |rV^ rcr))/V^. The two 
remaining eigenstates 

|T/i) = y^l^|0,0) + «i|5) (11) 

and 

1^2) = V^l-«||0,0)+t;2|5) (12) 
have eigenenergies 

En2) = l(s-2Ar± y/s^ + (5- 2A,)2^ (13) 

with 

+ ^i:ic-(^t.n)) 

r 

and 

^1(2) = = 

\jKh + - 2 A.) ((5 - 2 A, T v/8^e^ + - ^A.)^) 

(15) 

The existence of the K/ K' degree of freedom comphcates 
shghtly the definition of the CPS eigenstates. However, 
from the definition of \C±{t(j^ o^^^ can see that \S) 

corresponds to a generahzed spin-singlet state whereas 
I To), |T_) and |T+) correspond to generahzed spin-triplet 
states. The coupling teh hybridizes the empty state |0,0) 
with \S) only, due to the hypothesis that the supercon- 
ducting contact injects spin-singlet pairs inside the CPS. 
Figure 2. a shows the energies E2 and E_ as a func- 
tion of 5. The energies Ei and E2 show an anticrossing 
with a width 2\/2teh at J = 2A^, due to the coherent 
coupling between |0,0) and \S). The energy E- of the 
triplet states lies between Ei and E2. Figure 2.b shows 
the transition frequencies ujviT_^ ^t_V2 and UJV1V2 of the 
CPS, with ujm'm = {Em' —Em)/^. Thcsc frequencies will 
play an important role in the following. 

IV. MICROWAVE-INDUCED MATRIX 
ELEMENTS 

This section discusses the effect of the microwave gate 
bias on the eigenstates defined in section III. Inside the 



subspace 5, Hf^p has only three finite matrix elements, 
i.e. 

(T |Fff|Fi(2)) (16) 

and 

{T+ I i^^F 1^0) = ie— ^^(ofL^ac + (^R^ac) COs{ujRFt) 

(17) 

These terms are finite because Hf^p fiips the spins in the 
dots. The minus sign in Eq. (p!6|) is a direct consequence 
of the fact that |Vi(2)) comprises a singlet component 
whereas |T_) is a triplet state. In contrast, the plus 
sign in Eq. (pT|) is due to the fact that |To) and |T^) 
are both triplet states. The matrix element of Eq. (pT|) 
is always non-resonant since it couples two states with 
the same energy. Therefore, it can be disregarded in our 
study. The hamiltonian Hf^p has only one finite coupling 
element in the subspace f , i.e. 

{Vi I Hf^p \V2) = -viV2e{i^Lvt + i^RVac) M^Rpt) (18) 

with viV2 = y/2teh/\/^tlh + (^ - 2A^)2. The addition of 
f^LVac f^R^ac ^q. ([18]) is duc to the fact that the 
double occupation energy 6 is shifted by —{nhV^^it) + 
i^RV^(t)) when a microwave excitation is applied to the 
device. 

One can find experimental means to have V^^ and 
on phase, in agreement with the assumption made in sec- 
tion II. In this case, the matrix element (T_ | Hj^p |Vi(2)) 
vanishes when cxlv^^ = aRV^^. This effect is directly 
related to the injection of coherent singlet Cooper pairs 
inside the CPS since it is due to the existence of the 
minus sign in Eq. (p!6|) . If the injected pairs were in a 
product state instead of an entangled state, the matrix 
element (p!6|) would not be subradiant (see section VII. 4). 
Therefore, coherent pair injection inside the CPS can be 
revealed by observing microwave-induced transitions be- 
tween 1^1(2)) and and checking that these transi- 
tions are suppressed for aLV^^ = anv^^. The following 
sections describe how to probe these microwave-induced 
transitions with a DC current measurement. 



V. MASTER EQUATION DESCRIPTION OF THE 
CPS DYNAMICS 

In the following, the states |To) and |T+) are disre- 
garded because they are not populated in simple limits 
where relaxation towards them is neglected. The sequen- 
tial tunneling limit Fat <C ksT is furthermore assumed, 
with Fat the tunnel escape rate of an electron from one 
of the dots to the corresponding normal lead. For sim- 
phcity, it is assumed that this rate does not depend on 
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the dot orbital and spin indices. This would change only 
quantitatively the results shown in this paper. In the ab- 
sence of microwave irradiation, the dynamics of the CPS 
can be described with a master equation j30[ |4l[ 



— = MP 
dt 



with 



P = 



Pv. 
Pv, 
Pt_ 

-^single 



(19) 



(20) 



and 



M = 











(1 







-2w|rjv 





(1 


- vl)rN 








-2Tn 







2v^rN 


2v^Tn 


2rN 




-r^v 



(21) 



Above, Pi denotes the probability of state with 
i G {Vi,V2,T_}. The vector P also includes the global 
probability Pgingie of having a double dot singly occupied 
state. The use of this global probability is sufficient to 
describe the dynamics of the CPS because the single elec- 
tron tunnel rate F at to the normal leads is assumed to be 
independent from the dot orbital and spin indices. The 
various singly occupied eigenstates of H^qjj are defined 
in section VII. 3. The exact relation vf-\-V2 = l has been 
used to simplify the above expression of M. 

The microwave excitation Hf^^ can induce resonances 
between the states |Vi(2)) and |T_) while the excita- 
tion couples \Vi) and IV2). One can use a rotating 
frame approximation on independent resonances to de- 
scribe these effects. This approach is valid provided one 
of the microwave-induced resonance has a dominant ef- 
fect on the others, which requires the frequencies ojviT-^ 
ujt_V2 ^ViV2 to be sufficiently different. The rotat- 
ing frame approximation also requires to use small am- 

L(R) L(R) 

plitudes i^L{R)^ac and CiL{R)'^ac compared to ooviT-^ 
^T-V2i ^ViV2 and ojrf- In this case, the stationary state 
occupation probabilities can be obtained from 



= (M + MRF)Psi 



(22) 



with 



Mrf 



-rviT- - rviV2 

rViV2 

rviT- 




rviV2 
-rT^V2 - ^ViV2 
rT_V2 




rT_V2 
-rviT- - rT_V2 




(23) 






(24) 



CviT_ = 

Ct_V2 = 
CviV2 



Vie 
V2e 



2A 



(26) 

(27) 
(28) 



and ^^Pstat,i = 1- Above, Tab corresponds to the co- 
herence time between the states \a) and \b). Assuming 
that Tab is limited by tunneling to the normal leads, 
one obtains Tv.Vo = Pat, TviT_ = (1 + '^i)PAr and 



ViV2 

4) 



PT_y. ={1^v^)Tn. 

Figure 2.c represents schematically the dynamics of 
the CPS near the working point S = 2A^. Due to the 
assumptions made in section II on V5, the tunnel tran- 
sitions between the different CPS states (blue arrows) 
always occur together with the transfer of one electron 
towards one of the normal metal contacts. In contrast, 
the microwave irradiation induces transitions between 
the states |Vi) and II/2), |Vi) and |T_), or |T_) and {¥2) 
without any exchange of electrons with the normal con- 
tacts (red wavy arrows). The state |T_) can be reached 
through a microwave-induced transition but not through 
a tunnel process because it has no component in |0,0). 
The states |Vi) and IV2) can be both reached or left 
through a tunnel event because they have components 
in both |0,0) and \S). 

VI. RESULTS 
VI. 1 Principle of the measurement 

From Eq. (|2T]) . the tunnel rate transitions from the 
states IV2), and |T_) to the ensemble of the singly 
occupied states are 2viTn, 2^1 T at and 2FAr respectively, 
while the tunnel transition rate from a singly occupied 
state to or IF2) is Fat. As a result, the DC current 
/ flowing at the input of the CPS can be calculated as 



/ = R.P. 



stat 



(29) 



with R = eTN [2v'l , 2w| , 2, 1] . Fi sure 3. a shows the coeffi- 
cients v\ and ^2 as a function of 5. One can conclude from 
this plot that except dX 5 = 2A^, the various components 
of R have different values. Therefore, a microwave excita- 
tion changing the population of the states l^i), \V2)^ and 
|T_) should affect the value of the DC current flowing 
through the CPS. This effect will be used in the follow- 
ing to reveal the microwave-induced transitions between 
|F2>, and |T_). 



VI. 2 Stationary CPS state occupations 



rab{^) = 



\Cab 



2F 



ab 



> 



(25) 



Let us first discuss the dependence of the CPS state 
probabilities Pi on the parameter S for Vac^^ = (see 
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Fig. 3, black dotted lines in the three lowest panels). To 
understand this dependence, one must keep in mind the 
fact that the tunnel rate transitions from the states 
and to the ensemble of the singly occupied states are 
2vfTN and 2v2Tn^ as already discussed in section VI. 1. 
For 5 well below 2Ar, tends to zero. As a result, 
the CPS cannot escape easily from the state whose 
probability tends to 1. This is because in this limit, the 
state \Vi) is almost equal to the empty state |0,0), which 
makes the emission of an electron towards the normal 
leads very difficult. On opposite, for 5 well above 2A^, 
it is the probability of the state which tends to one 
because | V2) tends to |0, 0). In the absence of a microwave 
irradiation, the probability of state |T_) remains equal to 
zero since transitions towards these state are not possible. 
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FIG. 3: Coefficients '^1(2) 5 cind probabilities Pvi , Pv2 and Pt_ 
of the CPS states l^i), 1^2) and |T_) as a function of 6. We 



have used teh/^so = 1/3, Ak/k'/^so 



1.37 10 
3.10-^, 



3, 27r;iriv/A,o 



evi:JAso = 1/15, ev^JAso = 8/15, aL = an = 

K,L = K,R = 10~^ and UJRF — Steh- 



Let us now discuss the case Vac 

finite (see Fig. 3, 

red full lines in the three lowest panels). The term Hj^p 
excites the \Vi) ^ IV2) transition, which causes peaks or 
dips in Py^ and Pv2 for ujrf = i-^- S = Sy^^y^ 



2Ar ± 



(30) 



The term Hl^p excites the \Vi) ^ |T_) and |T_) ^ \V2) 
transitions, which causes peaks in Pt_ for uorf = ojviT- 
and ujRF = ^T-V2i i-^- ^ = ^Vi^T- and 5 = 5t_^V2 
respectively, with 



and 



5v,^T_ = 2Ar - UJRF + {2tl^/ujRF) (31) 



St_^V2 = 2Ar + OURF - {2tl^/uRF) (32) 



The term Hf^p also causes peaks or dips in Py^ and 
Py2, but they are hardly visible due to the scale used 
in Fig. 3. The decoherence rates Fy^T-, ^T-V2 and 
TviV2 have similar order of magnitudes (between Fat and 
2FAr). However, the width of the peaks or dips caused 
by Hj^p seems much larger than the width of the peaks 
caused by Hf^p. This is due to the limit (^l{r) ^ ^l{r) 
considered here. As long as the different types of res- 
onances are well separated in frequency, the resonance 
|Vi) ^ 1^2) gives probabilities Py^ and Pv2 which tend 
to the value 1/4 for ry^v2 sufficiently large. In princi- 
ple, the |Vi) ^ |T_) and |T_) ^ IV2) resonances give 
state probabilities Pi which saturate at more comphcated 
values which depend on 'u^2) when ry^T- and rT_V2 be- 
come sufficiently large. In the regime a^^R) <C i^l{r) 
considered here, the ^ 1^2) resonance is saturated 
while the \Vi) ^ |T_) and |T_) ^ IF2) resonances are 
only weakly excited. This explains that the width of the 
peaks or dips related to the \Vi) ^ IV2) resonance are 
much larger. 



VI. 3 Average current at the input of the CPS 

It is useful to discuss first the value /q of the current / 
at the input of the CPS in the absence of the microwave 
excitations. The current /q can be obtained from Eq. ([29|) 
with Vac'^^ = 0. From Fig. 4. a, /q shows a maximum 
for 5 = 2A^, where the two states and IV2) both 
correspond to equally weighted superpositions of |0,0) 
and \S). For S well below or well above 2A^, the current 
Iq vanishes because the CPS is blocked in the states \Vi) 
or IF2), respectively (see section V). 

Figure 4.b shows the difference between the current 
/ for a finite microwave irradiation and /q, as a func- 
tion of ujRF and S. The l^i) ^ 1^2) transitions yield 
a broad resonance along the curve urf = 0JV1V2 = 
V^^eh + ~ 2Ar)^ jh^ which has a frequency minimum 
ijjRF = 2\/2teh/^ at (5 — 2Ar. However, this resonance 
vanishes close to 5 = 2Ar because at this point, the 
tunnel escape rates 2v\Vn and 2'u|rAr of the CPS from 
\Vi) and IV2) are equal since vi = V2 and therefore, the 
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microwave- induced transitions between the states 
and IV2) cannot be seen anymore through a measure- 
ment of /. The |Vi) ^ |T_) and |T_) ^ \ V2) resonances 
yield two thinner resonances which cross at the point 
O corresponding to ^ = 2 A and Hujrf = V2teh' For 
ujRF tending to zero, the |Vi) ^ |T_) ( |T_) ^ \V2)) 
resonance progressively vanishes from / because this cor- 
responds to a regime where the state |Vi) (IV2)) is not 
populated anymore. Note that the calculation of the cur- 
rent / very close to the point O is in principle not valid 
using the rotating wave approximation on independent 
resonances since ojviT- = ^T-V2 this point. However, 
this represents only an extremely small area of Fig. 4. a 
(of order x Fat). Discussing the behavior of the CPS 
near point O goes beyond the scope of this paper. 



Figures. 4.c and 4.d show I — Iq as a function of S 
for two different values of ojrf- In Fig. 4.c, only the 
^ |T_) and |T_) ^ IV2) resonances are visible 
because ujrf < 2\^teh/^' In Fig. 4.d, the \Vi) ^ 
IF2) resonances are also visible. The ^ |T_) and 
|T_) ^ IF2) resonances appear as much thinner an 
smaller peaks. At the \Vi) ^ IV2) resonances, for the pa- 
rameters used in Fig.4.d, / reaches the saturation value 
Fat expected for rviV2 large and well separated reso- 
nances. This value can be obtained from Eq. (|29|) . using 

Pstat 

[1/4,1/4,0,1/2]. 



VI. 4 Dependence of the CPS input current on the 
amplitude of the microwave irradiation 

This section discusses how the minus sign in Equa- 
tion (p!6|) can be seen experimentally. One can note 
A/±^^^ = = ^vm) - ^0, A/v,^T_ = I{S = 
Svi^T-) - lo and AIt_^V2 ^T-^Vs) - ^0 the 

amplitudes of the microwave- induced current peaks ap- 
pearing for lurf = (^ViV2i ^RF = ^yiT_, and ujrf = 
UJT-V2' F)^^ to the symmetries of our model around 



the point 6 = 2A^, one has A/^ 



A/v 



Vi^V2 



AL 



Vi^V2 



and 



Lvi^T- = A/t_^v2- bottom panels of 

Fig. 5 show the variations of AIv^^T- and AIvi^V2 
with v^^ for a constant value of v^^. Due to the plus 
sign in Eq. ([18]), AIvi^V2 increases monotonically with 
v^^. In Fig. 5, this variation is very small because the 
^ 1^2) resonance is already saturated at v^^ = 
due to the value used for v^^. In contrast, due to the 
minus sign in Eq. ([T6]) , Aly^^T- shows a minimum for 
^ac — ^ac Note that if the electrons pairs injected in 
the CPS were not in an entangled state but in a prod- 
uct state, such a non-monotonic behavior would not be 
possible. For the parameters considered in Fig. 5, top 
panel, Aly^^T- vanishes at v^^ = v^^ because the effects 
of the \Vi) ^ IV2) resonance can be disregarded. This 



(b) 

(i-io)/er^ 




0.0 0.1 0.2 0.3 0.4 
I I I I I 



panel (d) 
^«RF=2f2t3, 




0.98 0.99 1.00 1.01 1.02 

5/2A 



0.90 0.95 1.00 1.05 1.10 

5/2A 



FIG. 4: (a): Current Iq in the absence of any microwave irra- 
diation as a function of S (b): Difference between the current 
/ for a finite microwave irradiation and the current /o, as a 
function of ujrf and 6 (c) and (d): Current difference I — Iq 
as a function of 5 for Hurf = 1.6te^ and Hurf = Steh- The 
other parameters used are the same as in Fig. 2. 



should not be true anymore in the case where the dif- 
ferent types of resonances are not well separated, which 
can happen e.g. if teh is too small with respect to the 
width of the resonances. However, in this case, one can 
still expect Aly^^T- to show a strongly non-monotonic 
behavior with a minimum at v^^ = v^^, provided the cou- 
plings aF{R) are sufficiently strong. Treating this case 
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requires to go beyond the rotating frame approximation 
with independent resonances used in this work. 

ev Va =2/3 



-z. 
U 

^ 0.01 


\ / 
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< 

0.1106 
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0.0 0.5 1.0 1.5 

V ^/v ' 



2.0 



FIG. 5: Amplitude of the current peaks AIvi^t_ (top panel) 
and AIvi^V2 (bottom panel) as a function of v^c for a con- 
stant value of Vac, i-e. eVac/^so = 2/3. The other parameters 
used are the same as in Fig. 2. 



VI. 5 Experimental parameters 

This section discusses the parameters used in the Figs, 
and the order of magnitude of the signals which can be 
expected in practice. In Figs. 3 to 5, the ratio of pa- 
rameters used correspond for instance to realistic val- 
ues teh = 50 /ieV (see Refs. ^so = 0.15 meV, 
^K/K' = 0.45 meV and Fat 



and Aly^^T- shown in Fig. 3.d are AIy^^v2 = 884 fA 
and Aly^^T- = 207 fA over a background Iq of 7.1 pA 
and 4.8 pA respectively. The maximum current difference 
A/y^^T_ in Fig. 5, top panel, corresponds to 129 fA for 
a background of 8 pA. Therefore, the features described 
in this article seem measurable experimentally. 



VII. DISCUSSION ON THE SPECTROSCOPIC 
ENTANGLEMENT DETECTION SCHEME 

The present section presents further examination and 
modifications of the model used above, in order to put 
the results of section VI into perspective. 



VII. 1 Use of a more general spin/orbit coupling 
term 

The Hj^p coupling term of Eq. ^ accounts for the cou- 
pling between the CPS and microwave excitations medi- 
ated by spin-orbit coupling. The above sections have 
used the particular form aira = icroii with G M, ob- 
tained from a microscopic description of spin-orbit cou- 
pling in a zigzag nanotube quantum dotjsij. This section 
discusses the generalization of the results to a more gen- 
eral coupling aira- Since Hf^p must be hermitian, one 
can use air^ = a^r and airi = o^*^ without any loss of 
generality. The parameter 

\Aa\ e'^" = {aLK + c^LK')Vac - i<^RK + (^rk')v^c (33) 

with ipo, g] — TT, tt] plays a crucial role in this case. It 
is convenient to redefine the states |T+) and |T_) more 
generally as 



with 



= 50 MHz (see Refs. [26|- and 
28|). Note that Fat = 50 MHz corresponds to 6.5 mK, 
therefore the sequential tunneling approximation used in I 
this work is valid using for instance T = 65 mK. In this 
case, using A = 5, the condition ([7]) to have electrons 
flowing only from the dot to the leads and not the re- 
verse gives V5 > 1.05 mV (see section II). This is com- 
patible with the condition 14 < A for having no quasi- 
particle transport between the superconducting lead and 
the dots, by using for instance a Nb contact for which 
A ~ 1.4 meV or a NbN contact for which A 3 meV. 
Using the above parameters, the ratio Vac^^ / ^so used in 
Figs. 3 and 4 corresponds to realistic microwave ampli- 
tudes v^^ — 10/iV and v^^ = 80/iV. Besides, the maxi- 
mum frequency urf = ^teh considered in this work (see 
Fig. 4.b) corresponds to 48.5 GHz, and the frequency at 
point O corresponds to 17 GHz^jwhich is accessible with 
current microwave technologies [42j. Using the above pa- 
rameters, the amplitude of the current peaks Aly^^y^ 



T+) = isgn((^,) (e-^^ 
T_) = isgn((^«) (e-^^ 



f^))/V2 (34) 
n))/V2 (35) 



-^{.^-l)\K.^K.) (36) 
-l(l + .^)|i.V,i.V) 

^^^{\Kcj,K'cj)^\K'cj,Kcj)) 

Note that |T+) and |T_) are still eigenstates of the hamil- 
tonian H^Jqj^^ with energy J— 2A^, corresponding to gen- 
eralized spin triplet states. The definitions of the other 
states 1^1(2)) and |To) remain unchanged. Using expres- 
sions (|34j) and (|35]) . one obtains 



(T+|if-^ 1^1(2)) = 



(37) 
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and 



For Ak^k' — ^ 0, one finds: 



= -'ievj |Aa| sgn{(pa) cos{LjRFt) 

for j e {1,2}. In sections II to VI, one uses 



(38) 



thus cpa = sgn(aL— Q^i?)7r/2 and \T±) 

wiiicii is in agreement witti Eqs.(|9]) and (p!Q|) . In ttiis limit, 
Eq. (|38l) agrees witii Eq. i^. Equations (|38l) and (|33l) 
siiow tiiat even witii a more general coupling term Hf^^^ 
the matrix elements (T_ | Hf^p |^i(2)) still present a sub- 
radiant form. Hence, the entanglement detection scheme 
discussed in this article appears to be quite general. Us- 
ing a more general Hf^^ will modify only quantitatively 
the predictions of section VI. 



VII.2 Role of Ak^k' / 



(T5I hso I Vi(2)) = le ^^___= ^ cos{ujRFt) (46) 



The coupling between | "1^1(2)) and IT5) is subradiant since 
it vanishes for o^rv^^ = o^L^ac- Nevertheless, for realis- 
tic parameters and in particular teh^ the transition 
frequencies cjTbVi andcj^bVi correspond approximately to 
2Ar/h^ which is too high for current microwave technol- 
ogy. This is why this paper focuses on microwave- induced 
transitions inside the subspace £. 



VII. 3 Microwave-induced transitions inside the 
singly occupied charge sector 

The different eigenstates of H^qjj in the singly occu- 
pied charge sector can be defined as: 



Remarkably, the subradiant matrix elements (p!6]) and 
(|38|) vanish for Ak^k' = 0. The aim of the present sec- 
tion is to show that using a finite Ak^k' does not rep- 
resent a fundamental constraint to have the subradiance 
effect. Indeed, |'i^i(2)) can still be coupled to other triplet 
states outside of the subspace £ when Ak^k' = 0. This 
fact is illustrated below, using a^r = io^i for simplicity. In 
this case, |Vi(2)) is coupled to a single triplet eigenstate 
\Tij) of H^Qjj outside the subspace defined by 











with 






_ ^K/K' 


~ 2Ar 


f2a) 


_ ^K/K' 


2Ar 



«- (l^it) - \T2i)) a+ (|r2t) - |Tu)) 



-al) 



(39) 



\K'a,K'a) - \Ka,Ka))+a-^ \K'a,Ka) 

(40) 

\K'a,K'a) - \Ka,Ka))+a^ \Ka,K'a) 



and 



a± = ctL - aR± —^{ceL + ceR 



Ar 



and 



such that 



One can check: 



cl{r) = o:L{Ry 



MR) 



K%D\Tb)=S\n) 



(41) 
(42) 

(43) 
(44) 



IM = yi -a^ (lifer, 0) - \0,Ka)) (47) 



2ArWl - a 



A, 



■.i\0,K'a)-\K'a,0)) 



|ai.) = -IJl - {\Ka,0) + \0,Ka)) (48) 



{\0,K'a) + \K'a,0)) 



2Ar\ l-(J- 



\h2.) = - W 1 + d^'^' 0) - |0' Ka)) (49) 



2A^^/l + (j^ 



{\Q,K'cj)}-\K'a,Q)) 



\a2.) = \Jl + cj^{\Kcj,0) + %Ka)) (50) 



+ 



2ArJ\ + (J^ 



(|o,ii:V)} + |ifV,o)) 



for (7 G {t^i}- These states have eigenenergies ej, 
and £2 respectively, with 



and 



(^6 I ^so \y\{2)) = ^eV^2)-^{^R - ^l) COs{ujRFt) (45) 



e',=e-tee^{-iyAr 



=e^tee^{-iyAr 



(51) 
(52) 



At 



-{aL - ^rY 



for i G {1,2}. The states \hia) and \h2c7) can be seen as 
generalized bonding states and |aicr) and |a2a) as gener- 
alized antibonding states. This section uses air — 
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for simplicity. The term couples |6^cr) and \aicr) 

to \bia) and \aia) only, for i G {1,2}. Only the transi- 
tions \aicr) ^ Ibia^) correspond to a finite frequency, i.e. 
ujai^bi^ = 2tee/^' One can check that 

whereas 



pied state with a spin cr{a^) on dot L{R). In the case of 
coherent Cooper pair injection, the double quantum dot 
effective hamiltonian can be written 1301 



Hdqd =£{nLa^nR^) 



ie{aLV^^-aRV^J cos{uRFt)/2 (53) 



(55) 



+ (te JV^) (4^4^ - 4^4^ + /i.e.) 



\HI?p\h 



\hia) 



\aia) 



(54) 



H 



RF 



Importantly, the matrix element of Eq. (|53|) has a sub- 
radiant structure. This property is due to the fact that 
the states and \aicr) are entangled states with differ- 
ent symmetries, i.e. \bi-^) is an antibonding state which 
contains some |ra, 0) + |0, ra) components whereas \aia) 
is a bonding state which contains |rcr, 0) — |0,rcr) com- 
ponents. This is analogous to the fact that the elements 
(T_ I Hj^p \ Vj) couple a state \ Vj) with a spin-singlet com- 
ponent to a spin-triplet state In contrast, the ma- 
trix elements of Eq. (|5^ are not subradiant because they 
couple two entangled states with similar symmetries, i.e. 
two bonding or two antibonding states. 

Due to the subradiant form of Eq. (j53j), the transi- 
tions \ai(j) ^ \bia) can lead to a non-monotonic varia- 
tion of the CPS input current as a function of e.g. v^^, 
due to another type of entanglement than the one dis- 
cussed in section VI. Therefore, in the context of the 
characterization of split Cooper pairs entanglement, one 
needs to find a way to discriminate possible current res- 
onances corresponding to the transitions \ai(j) ^ \bi-^) 
and 1^1(2)) ^ (^-|- practice, this should be feasi- sis A = {Vi,V2, 
ble by studying how the different resonance frequencies 
vary with the DC gate voltages of the two dots. Indeed, 
uJai^bi^ does not depend on the parameter ^, contrarily to 
^ViT- and ujt_V2- Therefore, possible current resonances 
due to \ai(j) ^ \bi-^) transitions should appear as hori- 
zontal lines in Fig. 4.b. This effect was nevertheless dis- 
regarded in section VI, assuming that ujai^bi^ is too large 
to be accessible experimentally. Studying quantitatively 
the possibility to observe the resonances \aia) ^ \bia) 
requires to go beyond the approximation of an electronic 
tunnel rate F at to the normal leads which is independent 
from the dot orbital and spin indices [43| . 



where Hint still forbids the double occupation of each 
dot. One uses above riia- = d\^dia- with the creation 
operator for an electron with spin a in dot i G {L^R}. 
Let us furthermore assume that there also exists a spin- 
flip coupling term to the microwave signal, with the form 



- dl^di^) 



(56) 



Such a spin-flip coupling can be due for instance to the 
magnetic field associated with the microwave irradiation. 
In practice, this term should have a very weak amplitude, 
but it is nevertheless discussed here for fundamental pur- 
poses. 

The term in teh hybridizes the CPS empty state |0,0) 
with the singlet state 5^ = (|t,|) - ||,t))/V^, so 
that an anticrossing appears again in the spectrum of 
the CPS even-charged states. For simplicity, it is ass- 
sumed below that the double occupation energy S = 2e 
of the CPS is degenerate with the energy of |0,0), i.e. 
^ = 0. In this case, one can use the orthonormalized ba- 



T 



Tb 



To)} of eigenstates of (j55j) 



in the even charge sector , with 



(2) 



(|0,0) 



± 



s 



and 



Ta(6)) = (|t,t)±|;,i))/V2 

ro} = (|t,i> + i;,t»/V2 



(57) 
(58) 

(59) 



The states 



Vi ) and 



V2 ) have energies Ei and E2 given 



by £^1(2) = ^teh' They play the role of the states \Vi) 



VII.4 Simplified model without the K/K' degeneracy and {¥2) of section VI. It is convenient to define 



Ta) and 



It is interesting to discuss a model without the K/K' 
degree of freedom to show simply how the subradiance 
property arises. 

Case of coherent Cooper pair injection 

Let us assume that each of the two CPS dots has a 
single orbital. One can note |cr, cr') a CPS doubly occu- 



rs ) as superpositions of triplet states with equal spins. 



The states 



To 



Ta) and 



T5) have an energy S 



0. 



One can check straigthforwardly that 



H 



sf 
RF 



S) = {Xr-Xl) 



thus 



n 



H 



RF 



1(2) 



±{\r - Al)/2 



(60) 
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whereas { T. 



a(0) 



rrsf 



(2) 



and {V2 



H 

sf 



sf 
RF 



0. 



The states Vi(2) are thus coupled by H^p to a single 
state 75^, with a subradiant matrix element (|6Q|) . This 

illustrates the universality of the mechanism discussed in 
section VI. 



are expected for a model including the K/K' degree of 
freedom. This illustrates that the subradiance property is 
a good indication of the injection of entangled Cooper 
pairs inside the CPS. More sophisticated descriptions of 
incoherent injection of Cooper pairs into the CPS are 
beyond the scope of this article. 



Case of incoherent singlet injection 

One can model naively the incoherent injection of 
Cooper pairs inside the CPS by assuming that up spins 
are always injected inside the left dot and right spins 
inside the right dot. This requires to replace the hamil- 
tonian ([55]) by 



DQD = ^{riLa ^ riR^) 

+ teh{d]^^d\^^ + dR^dL^) + Hint 



(61) 



One can use again (5 = 26: = for simplicity. In this 
case, one can define an orthonormalized basis B = 



{Wi,W2, 



Tg ) } of eigenstates of (f6T]) in the 



even charge sector, with 



|M/i(2)) = (|0,0)±|t4))/V2 



(Afi |t,t) + AL |;,;))/\/Ai + A^ 



(AL|t,t)-Afl \i,\-))I^L 



(62) 



(63) 



(64) 



and 



Te^ = |;,t). The role of the states \Vi) and \V2) of 
section VI is now played by \Wi) and IW2). The states 
\Wi) and 11^2) have again energies Ei and E2 defined in 



the previous section, whereas the states 
Tp ) have an energy (5 = 0. The states 



Td) and 



and 



of the previous section are still CPS eigenstates, but it 

fr) and 



is more convenient to use the eigenstates 

to study the effect of H^^. Due to the term in tg/i, the 
states \ Wi) and \ W2) still form an anticrossing in the en- 
ergy spectrum of the CPS. Hence, such an anticrossing is 
not characteristic from the injection of entangled Cooper 
pairs. The only state of B connected to |Wi(2)) by 

is Tc ) , with a matrix element 



VIII. COMPARISON WITH AN ALTERNATIVE 
SETUP: THE CPS EMBEDDED IN A 
MICROWAVE CAVITY 

This section compares the experimental scheme pro- 
posed in Ref.[22j to the scheme discussed in the present 
article. Reference [22[ suggests to observe the minus sign 
in Equation ([16]) by inserting the CPS inside a coplanar 
microwave cavity to obtain a lasing effect involving the 
|Vi) |T_) transition. The minus sign in Equation ([T6|) 
leads to a non-monotonic dependence of the number of 
photons in the cavity as a function of the coefficients aL 
and au^ which mediate a coupling between the CPS and 
the electric field conveyed by the cavity. Since this elec- 
tric field can be considered as constant over the whole 
CPS area, it is necessary to be able to vary aL indepen- 
dently from aR to observe a non-monotonic behavior in 
the number of photons. This can require to complexify 
the CPS design, for instance. In the present scheme, such 
a control on aL and aR is not necessary since it is suffi- 
cient to vary independently the amplitudes v^^ and v^^. 
This can be naturally achieved by using two independent 
microwave supplies for the two gates. The advantage of 
the scheme presented in Ref.[22] is that the signal to be 
measured is a large photon number which can be ob- 
tained by measuring the power spectrum emitted by the 
cavity. In other words, the scheme of Ref.[22] exploits 
the fact that the lasing effect provides an intrinsic am- 
plification process for the Vi ^ T_ transitions. In the 
present scheme, the measurement seems a bit more dif- 
ficult since the current peaks to be measured are very 
small. Nevertheless, such current amplitudes are mea- 
surable, in principle [42]. Therefore, the scheme presented 
in this reference could be an interesting alternative ap- 
proach to demonstrate the coherent injection of singlet 
Cooper pairs inside a CPS. This scheme furthermore al- 
lows to study also the {¥2) ^ (T_ | transition, which is 
not possible with the scheme of Ref. |22|] . 



IX. CONCLUSION 



H'jF\Wm) = ±x/2{Xl^X%) 



(65) 



which is not subradiant, but increases monotonically with 
Xr and Xl- Therefore, the subradiance property is lost 
when Cooper pairs are injected inside the CPS in a prod- 
uct state instead of an entangled state. Similar results 



The DC current response of a double-quantum-dot 
based Cooper pair beam splitter (CPS) to a microwave 
gate irradiation is a very rich source of information on 
Cooper pair splitting. In particular, it can reveal the 
entanglement of spin-singlet Cooper pairs injected inside 
the CPS. This article illustrates this property for a double 
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quantum dot formed inside a carbon nanotube with typ- 
ical parameters. If they are spin-entangled, the injected 
pairs are coupled to other CPS states through some 
subradiant microwave transitions mediated by spin-orbit 
coupling. This property can be revealed by applying to 
the two CPS quantum dots two on-phase microwave gate 
voltages. The spin-orbit mediated microwave transitions 
cause DC current resonances at the input of the CPS. 
The subradiance property manifests in a strongly non- 
monotonic variation of these current resonances with the 
amplitude of the microwave signal applied to one of the 
two CPS dots. This behavior does not depend on details 
of the model like the exact form of the spin-orbit interac- 
tion term. Similarly, the presence of atomic disorder in 
the nanotube has to be assumed only for quantitative rea- 
sons. More generally, the entanglement detection scheme 
discussed in this work could be generalized to other types 
of quantum dots with spin-orbit coupling like e.g. InAs 
quantum dots, in principle. For simplicity, this article 
discusses the limit where the intra-dot charging energies 
are very strong, so that there cannot be two electrons 
at the same time on the same dot. For smaller charging 
energies, the efficiency of Cooper pair splitting should be 
decreased. Nevertheless, if the CPS produces entangled 
split Cooper pairs with a sufficient rate, the resulting sub- 
radiant current peaks should still be observable. Inter- 
estingly, the bonding or antibonding single particle states 
delocalized on the two dots of the CPS can also cause a 
subradiant current resonance, because they present an- 
other type of entanglement. However, in principle, this 
resonance can easily be discriminated from the subradi- 
ant resonances caused by split Cooper pairs, because of 
a different dependence on the CPS DC gate voltages. 

/ acknowledge useful discussions with T. Kontos, J. 
Viennot and A. Levy Yeyati. This work has been financed 
by the EU-FP7 project SE2ND[271554]. 
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